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We address the use of entangled qubits as quantum probes to characterize the dynamical noise induced by 
complex environments. In particular, we show that entangled probes improve estimation of the correlation time 
for a broad class of environmental noises compared to any sequential strategy involving single qubit preparation. 
The effect is present when the noise is faster than a threshold value, a regime which may always be achieved 
by tuning the coupling between the quantum probe and the environment inducing the noise. Our scheme ex¬ 
ploits time-dependent sensitivity of quantum systems to decoherence and does not require dynamical control 
on the probes. We derive the optimal interaction time and the optimal probe preparation, showing that it corre¬ 
sponds to multiqubit GHZ states when entanglement is useful. We also show robustness of the scheme against 
depolarization or dephasing of the probe, and discuss simple measurements approaching optimal precision. 
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The coherence properties of a quantum system are strongly 
affected by its interaction with the surrounding environment. 
This is often an obstacle to the implementation of quantum 
technologies, so that much effort has been devoted to study 
the system-environment interaction and to engineer decoher¬ 
ence in order to minimize its degrading effects l[I]|2l. On the 
other hand, the very sensitivity of quantum systems to exter¬ 
nal influence also provides an effective tool to characterize 
unknown parameters of a given environment Ea by exploit¬ 
ing quantum probes, as opposed to classical ones, usually big¬ 
ger and more intrusive. Indeed, characterizing the dynamical 
noise induced by an external complex system is of great rele¬ 
vance in many areas of nanotechnology, as well as in mon¬ 
itoring biological or chemical processes ESI. Besides, it 
represents a crucial step to design robust quantum protocols 
resilient to noise EHH. 

The proper framework to address characterization by quan¬ 
tum probes iniiia, and to design the best working condi¬ 
tions, is given by quantum estimation theory iflTlI . which pro¬ 
vides analytical tools to optimize the three building blocks of 
an estimation strategy: (i) preparation of the probe system in 
a suitably optimized state, (ii) controlled interaction of the 
probe with the system for an optimal amount of time t, (iii) 
measurement of an optimal observable on the probe. Over¬ 
all, the ultimate precision for any unbiased estimator 7 of a 
certain parameter 7 is bounded by the quantum Cramer-Rao 
(CR) theorem, stating that Var( 7 ) > [MiT( 7 )]“^, where M is 
the number of measurements and H ( 7 ) is the quantum Fisher 
information (QFI), i.e. the superior of the Fisher information 
over all possible quantum measurements described by positive 
operator-valued measures (POVMs). 

Recently, single-qubit quantum probes have been proposed 
for the characterization of dynamical noise by monitoring de¬ 
coherence and dephasing induced by the environment under 
investigation lfT8l - E2l . In this paper, we extend this analysis 
to entangled qubits used as quantum probes, and show how 
they greatly improve the characterization of a broad class of 
environmental noises compared to any sequential strategy in¬ 
volving single qubit preparation Il23l425l . In particular, we 


show how to improve estimation of the correlation time (i.e. 
the spectral width) of classical dynamical noise described a 
stochastic process. Indeed, stochastic modeling of the envi¬ 
ronment has been proved an effective tool to reliably describe 
dynamical noise and the decoherence process in several sys¬ 
tems of interest Il26l434l . including motional averaging iTSll 
and solid state qubits Il36l438ll . In particular, since such noise 
is usually emerging from a large collection of fluctuators, we 
are going to consider Gaussian stochastic processes. 
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FIG. 1. Schematic diagram of possible characterization techniques. 
On the left: we have N qubits and each one is prepared in the state 
|-|-), then interacts with the system for a time t and is finally mea¬ 
sured independently of the other qubits. On the right: the N qubits 
are initially prepared in a multiqubit GHZ state |T'n), and are let 
interact with the system such that they are subject to the same re¬ 
alization of the noise. At the output, a collective measurement is 
performed on the qubits. 

The probing scheme is depicted in Fig. [T] both for a se¬ 
quence of N uncorrelated qubits and for an W-qubit entan¬ 
gled state. In both cases we assume that the qubits do not 
interact with each other. We also assume that in the first 
case each qubit is subject to different realizations of the noise, 
while all the qubits in the entangled state interact with a com¬ 
mon environment. Let us start by considering a single qubit 
interacting with a fluctuating dephasing environment. The 
Hamiltonian is given by T-L{t) = -f B{t)(Jz, where 

Wo is the energy of the qubit and B{t) is a realization of 
the stochastic process describing the dynamical noise. As 
a paradigmatic example we consider a zero-mean Ornstein- 
Uhlenbeck process characterized by the autocorrelation func¬ 
tion Ar(f,T) = {B{t)B(t'))g = i 7 rexp[— 7 |f — f'|], or by 
the corresponding Lorentzian spectrum. Here, 7 is the spec- 
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tral width, i.e. the inverse of the autocorrelation time, while F 
denotes the coupling between the probe and the system. It is 
worth noticing that a similar analysis may be carried out for 
other Gaussian processes, e.g. processes with power-law or 
Gaussian autocorrelation functions, and that results are quali¬ 
tatively the same, independently on the choice of the autocor¬ 
relation function. 

The density operator of the evolved qubit is given by 

pit) = {Uit)piO)U\t))s, ( 1 ) 

where {■)g denotes the average over all possible realizations 
of the stochastic process in the time interval [0,f], U{t) = 
exp[—i Jg 'H{s)ds] = eyip{[—iuJot + 4’{t)\cjz} is the time evo¬ 
lution operator, and (j){t) = J* B{s)ds is the accumulated 
phase due to the environmental (dynamical) noise. An explicit 
expression for p{t) can be found by employing the charac¬ 
teristic function of a zero-mean Gaussian stochastic process; 
did) where 

= [ [ ds dw K{s,w) =+'jt — 1 ). (2) 
Jo Jo 1 

If the qubit is initially prepared in a state described by the 
density operator p(0), the density operator at the time t will 
be p{t) with pkk{t) = pfefc(O), fc = 1, 2 and 

Pi2(i)=e-2[-o‘+/3.(*)lpi2(0). (3) 

The optimal single qubit preparation, given by |-|-) = 

^(|0) -f |1), together with the cotTesponding QFI and the op¬ 
timal measurement for the estimation of 7 have recently been 
found II 22 I . For N uncorrelated qubits, thanks to additivity, 
the QFI is just N times the single qubit QFI, i.e. 

= ( 4 ) 

Let us now consider a probe made of N qubits initially pre¬ 
pared in the generalized GHZ entangled state 

= + (5) 

We also assume that the qubits do not interact with each other 
and that they interact with a common environment, i.e. with 
the same realization of the noise. The overall Hamiltonian is 
thus 

•HW(f) = n{t) (g) -f I (g) n{t) (g -h ..., (6) 

where 'H(f) is the above single qubit Hamiltonian and B{t) is 
the same realization of the noise for all the qubits. The QFI 
for the parameter 7 reads 

_ I [^7^7(0] , ( 7 ) 

where P-y {t) (see Eq. (|^) is a monotonically increasing func¬ 
tion of t with Pj{0) = 0. Moreover, we have ^ 


r(f—1/7) forf 1 . Thus both (7, f) andiF™^(7,t) are 
asymptotically vanishing and show a single maximum, corre¬ 
sponding to different optimal values of the interaction time, 
and respectively. We refer to this maximum value 
as the maximal QFI for a specihc value of 7. As is appar¬ 
ent from the above equations, the maximization of the QFI 
involves transcendental equations and must be done numeri¬ 
cally. The behavior of t) and t) is depicted in 
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FIG. 2. (Color online) The QFI ratio Rn as a function of 7 (log scale) 
for A = 2 (solid red), A = 3 (dashed blue) and A = 4 (dotted 
green). For small values of 7 , the ratio is below one (black dashed 
line), and tends asymptotically to 1/A: in this regime it is more con¬ 
venient to employ separable states than maximally entangled states. 
The ratio increases monotonically with 7 and exceeds one at a thresh¬ 
old value 70 (A), which depends on A. For asymptotically large 7 , 
the use of A-qubit GHZ states is A times better than the use of A 
qubits in a separable state. In the inset we show the QFI 
(solid) and (dashed) for 7 = 10 . 


In situations where we can control the interaction time be¬ 
tween the probe and the environment, it will be most conve¬ 
nient to set it to the optimal time. Thus, a fair comparison be¬ 
tween separable and entangled probes naturally involves the 
maximal QFI of the two cases. We therefore introduce the 
QFI ratio as 


Rn(7) 




(8) 


and analyze its behavior as a function of 7 and N. When 
Rn( 7 ) > 1 , the use of a A qubits in a GHZ state improves es¬ 
timation compared to the use of N uncorrelated probes, e.g. in 
a sequential strategy. Figure [^illustrates the main results: the 
ratio Rn( 7) is larger than one for 7 > 70(A), where 70(A) is 
a threshold value that depends on A. Moreover, Rn( 7) —f A 
for 7 70(A). This result is enhanced by the fact that, upon 

substituting 7 = 7/r and t = pt, we may show that the 
quantum signal-to-noise ratio (QSNR) does not de¬ 

pend on F. This means that, if one is able to control the cou¬ 
pling between the probe and the system, one can always tune 
F to achieve a situation where Rn( 7) > L All the Figures are 
obtained by setting F = 1 . 
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Now a question arises: Is the GHZ state |^'n) the optimal 
one? Are there (entangled) states that give even higher QFI? 
The answer to this question cannot be analytic, because one 
cannot diagonalize analytically a generic density matrix of a 
multiqubit state. In order to attack this problem, we first no¬ 
tice that the maximum is achieved for an initial pure state ll24l . 
We have thus generated a large number (n = 10 ®) of random 
initial pure states, uniformly distributed according to the Haar 
measure, for different values of 7 and for N = 2,3,4. For 
each random state, the maximal QFI, resulting from 

the interaction with a common environment has been numeri¬ 
cally evaluated using the expression 


H{-,)=2J2 


(V' m Id^Pjllpn) 


(9) 


where |?/;„) (-(/inl is the diagonal form of the den¬ 

sity operator after the interaction with the environment. This 
value is then used to evaluate the corresponding QFI ratio 
{'j) /{'-f), and to compare the estimation precision to 
the precision achievable using N independent qubits interact¬ 
ing with separate environments. Our results indicate that, for 
7 ^ 7 o^ that is, in the region where entanglement is conve¬ 
nient, the GHZ state is indeed the optimal one, thus showing 
that entanglement is a resource for the estimation of the spec¬ 
tral width of Gaussian dynamical noise. Below the threshold 
the GHZ state is no longer optimal and there is numerical ev¬ 
idence that an extremal state can be found by optimization in 
the same family of states that had been identified in |[39l as 
optimal probes to improve frequency estimation (notice, how¬ 
ever, that below threshold 7 ^ 70 , the optimal strategy is the 
use of uncorrelated qubit probes interacting with independent 
environments). The states of this family, for N qubits, have 

the form 1 $^) = ^k\k), where are normalized real 

coefficients, \k) is an equally weighted superposition of all N- 
qubit states with a number fc or a number N — k of excitations, 
and [ J denotes the integer part. The GHZ state belongs to this 
family with oq = l/v /2 and all other coefficients set to 0 . 

Figure illustrates our numerical results obtained for two 
and three qubits. The plots show the QFI ratios R 2 (left) and 
R 3 (right). The solid blue line is the ratio for the GHZ state, 
the 10^ gray points correspond to the QFI ratio of randomly 
generated states and the dashed red line is found by optimizing 
the QFI over the coefficients Ofc of |4 )n). We can see that from 
7 ~ lo{N) the blue and red curves coincide, i.e. GHZ states 
are extremal. We also notice that for 7 > 7 o(^) a significant 
fraction of gray dots lies above the r = 1 dashed line, but the 
dots are sparse around the solid blue line, meaning that the 
GHZ state allows for a remarkable gain in the estimation of 
larger values of 7 . 

It is worth to emphasize that optimal precision, i.e. the QFI 
of Eq. (|7]| may be achieved upon implementing a simple rank- 
2 measurement. Indeed, corresponds to the Fisher 

information of a projective measurement on the two eigenvec¬ 
tors corresponding to the nonzero eigenvalues of the evolved 
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FIG. 3. (Color online) The curves in the logplots show the QFI ra¬ 
tio R2 (top) and R3 (bottom) as a function of 7 for the GHZ state 
(solid blue) and the optimal state of N qubits in a common environ¬ 
ment (dashed red). The two curves superimpose each other above 
the threshold 70 (N). We also show the QFI ratio for 10° randomly 
generated states (gray points), uniformly distributed according to the 
Haar measure. 


density operator, which are, respectively 

i (1 ± (10) 

|p±) = ± _ (11) 

The CR theorem, however, sets only a lower bound to the pre¬ 
cision of any unbiased estimator, and a question arises on how 
to suitably process data coming from the above rank -2 mea¬ 
surement in order to saturate the bound. Bayesian estimators 
are known to saturate the CR bound for asymptotically large 
numbers of measurements: in order to assess quantitatively 
the performance of Bayesian estimation we have performed 
simulated experiments on the probing system. In particular, 
we have simulated the outcomes {xi ,..., xm} of the mea¬ 
surement by randomly choosing a result according to the prob¬ 
abilities of Eq. ( [Tol l have built a Bayesian estimator 7 as 
the mean value of the a posteriori distribution, starting from 
a flat prior. The resulting relative error of 7 , e = ■^Var( 7 )/ 7 , 
is shown, as a function of the number of measurements, in 
Eig. for a specific value of 7 . We see that with a relatively 
low number of measurements, of the order of thousands, the 
bound is saturated and the situation improves by increasing 
the number of qubits. The proposed scheme thus allows for 
an effective and achievable estimation of the parameter 7 . 

Let us now address the robustness of our scheme against 
noise in the preparation of the probe. In fact, we have shown 
that the use of entangled qubit probes prepared in a GHZ state 
leads to enhanced precision in the estimation of the spec¬ 
tral width. However, it is generally challenging to exper¬ 
imentally prepare the probes exactly in the GHZ state and 
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FIG. 4. (Color online) Log-log plot of the relative error e = 
^Var( 7)/7 of the Bayesian estimator as a function of the number 
of measurements, for 7 = 10. The lines represent the CR bound 
for a single-qubit measurement (dashed blue) and for a 4-qubit GHZ 
state (solid red). The diamonds (single qubit) and the dots (4 qubits) 
correspond to the performance of a Bayes estimator applied to sim¬ 
ulated experiments. Bayes estimators saturate the CR bound when 
increasing the number of measurements and very good performances 
are achieved already for thousands of measurements. The plot also 
shows that estimation improves with the number of qubits since the 
CR bound is saturated with a lower number of measurement. 


a question arises on how sensible is this estimation scheme 
to, e.g., the purity of the initial preparation. We answer 
this question by considering a partially depolarized state, 
Pp — P Pghz + (1 ~ where I is the identity ma¬ 

trix and 0 < p < 1, and a partially dephased state, pg = 
(5poHz+5(l-<5)(|0)®^(0|-|-|l)®^(l|)where0 < (5 < 1. In 
both cases, an analytic expression for the QFI may be found; 
we have 




2tv+2jv4 [(2^_2 )p + 2] p^[d,l3^{t)]^ 
i(2^ - 2) p + 2f e4^"/3(t,7) _ 4Np2 

g47V2/3(t,7) _ ^2 [^ 7/^7 (^)] > (13) 


which are obviously less than being pp and pg mixed 

states, but may be still larger than Indeed, Figure 

shows that for each value of 7 above ‘jo{N) there is a thresh¬ 
old value for the purity, above which the use of a depolar¬ 
ized or dephased GHZ state still leads to an improvement over 
the use of N uncorrelated probes. The threshold purity po is 
close to one for 7 ~ 'yo{N) and for 7 ^ 7o(A^), whereas it 
shows a minimum in the intermediate region, thus allowing 
for a certain tolerance in the preparation of the initial state of 
the probe. Besides, this minimum value of the threshold gets 
lower when increasing the number of qubits. 

In conclusion, we have shown that the use of entangled 
qubits as quantum probes outperforms the sequential use of 
single-qubit probes in the characterization of the dynamical 
noise induced by complex environments. In particular, we 
have shown that entangled probes improve estimation of the 
correlation time for a broad class of environmental noises 
when the noise is faster than a threshold value. This result 



FIG. 5. Threshold purity po, as a function of the spectral width 7 , 
above which the use of a depolarized (left panel) or a dephased GHZ 
state (right panel) is more effective than a set of uncorrelated qubits. 
The different lines correspond to a different number N of qubits: 
N = 2 (solid), N = 3 (dashed), = 4 (dotted), N = 5 (dot- 
dashed). The threshold po approaches 1 when 7 —>■ 70 (A^) and 
when 7 —^ 00 , whereas there is an intermediate region where po 
decreases to a minimum, meaning that there is more tolerance in the 
initial preparation of the probe. When N increases the minimum of 
Po decreases and moves to larger values of 7 . 


is enhanced by the fact that, upon controlling the coupling be¬ 
tween the probe and the system, the threshold value can be 
reduced arbitrarily. Our scheme exploits time-dependent sen¬ 
sitivity of quantum systems to decoherence and does not re¬ 
quire dynamical control on the probes. We have found the op¬ 
timal interaction time and the optimal multiqubit probe prepa¬ 
ration, showing that it corresponds to multiqubit GHZ states. 
The proposed measurement scheme achieves the Cramer-Rao 
bound for a relatively low number of measurements, upon em¬ 
ploying a Bayesian estimator, and is robust against imperfect 
preparation of the initial entangled state. 
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